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Understanding the metal-insulator transition in disordered many-fermion systems, both with and 
without interactions, is one of the most challenging and consequential problems in condensed matter 
physics. In this paper we address this issue from the perspective of the modern theory of the 
insulating state (MTIS), which has already proven to be effective for band and Mott insulators in 
clean systems. First we consider noninteracting systems with different types of aperiodic external 
potentials: uncorrelated disorder (one-dimensional Anderson model), deterministic disorder (Aubry- 
Andre Hamiltonian and its modification including next-nearest neighbour hopping), and disorder 
with long-range correlatioirs (self-affine potential). We show how the many-body localisation teirsor 
defined within the MTIS may be used as a powerful probe to discriminate the insulating and 
the metallic phases, and to locate the transition point. Then we investigate the effect of weak 
repulsive interactions in the Aubry-Andre Hamiltonian, a model which describes a recent cold- 
atoms experiment. By treating the weak interactions within a mean-field approximation we obtain 
a linear shift of the transition point towards stronger disorder, providing evidence for delocalisation 
induced by interactions. 

PACS numbers: 71.10.Fd,71.30.+h,71.23.An,72.15.Rn 


I. INTRODUCTION 

In the modern theory of the insulating state (MTIS), 
which was initiated by the seminal article written by 
Kohiii, the different behaviours of metals and insula¬ 
tors are attributed to the different organisations of the 
electrons in the many-body ground state. In insulators, 
the electrons satisfy a many-particle localisation condi¬ 
tion^. Kohn associated this localisation with the break¬ 
up of the many-electron wave-function into terms which 
are localized in essentially disconnected regions of the 
many-particle configuration space. This approach is fun¬ 
damentally different from conventional theories of insula¬ 
tors, which require the knowledge of, at least, low-lying 
excitations, and are tailored towards a specific kind of 
insulator, depending on the physical mechanism which 
triggers the insulating behaviour. Some fundamental de¬ 
velopments in the MTIS were achieved only in the late 
1990’s (thanks to the works by Resta and Sorellai^, and 
others^i^). These developments began with the observa¬ 
tion that the polarization is finite in insulator, while it 
is ill-defined in metals. This led to a definition of the 
localisation tensor— of a many-body system which is de¬ 
rived from the Berry-phase formulation of the polarisa¬ 
tion as understood within the modern theory of the po- 
larization^^— . The direct connection between the many- 
body localisation tensor and the disconnected parts of 
the many-electron wave-function — as originally defined 
by Kohn — was demonstrated by Souza, Wilkens and 
Martins^. 

The MTIS is supposed to be adequate to describe any 
kind of insulator, independent of the physical mechanism 
which induces the insulating behaviour. It should apply 
to band, Mott^, Anderson^, quantum-Halli^ and possi¬ 
bly even to topological insulators. Therefore it repre¬ 


sents a promising approach to address the outstanding 
open problem of the fate of Anderson localisationii in 
the presence of interactionsi^rJi. 

So far, band and Mott insulators have been analysed 
within the framework of the MTIS^, both using lattice 
models in a tight-binding scheme, and also via ab initio 
electronic structure simulations^d^^J^. Instead, Ander¬ 
son insulators have received little attention. In particu¬ 
lar, it is not known — even for the noninteracting case — 
whether the many-body localisation tensor may be used 
to locate the critical point of the (Anderson) transition 
which separates the conducting and the insulating phases 
in disordered systems. 

The first purpose of the paper is to investigate this issue 
in the noninteracting case. With this aim, we study the 
Anderson transition in noninteracting one-dimensional 
lattice models close to half filling. Since disorder cor¬ 
relations play a fundamental role in low dimensional sys¬ 
tems (for instance, they determine the presence or ab¬ 
sence of transition points and mobility edgesiSr^^), we 
consider various models of disorder and study the many- 
body localisation tensor, and its reliability, in qualita¬ 
tively different scenarios. First, like Ref. |^, we con¬ 
sider one-dimensional lattices with uncorrelated disorder, 
where the single-electron orbitals are expected to be ex¬ 
ponentially localised at any nonzero disorder strengtbii. 
Then we focus on the more intriguing and instructive 
case of deterministic disorder due to an external peri¬ 
odic potential whose period is incommensurate with the 
lattice. All single-particle orbitals of this Hamiltonian 
(named Aubry-Andre model) become localised, but only 
beyond a finite disorder strength^^. Next we consider a 
generalised Aubry-Andre model including next-nearest 
neighbour hopping, where a mobility edge separating 
extended and localized single-particle orbitals was pre- 
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dicted^^. Further, the case of non-deterministic disorder 
with tunable spatial correlations is addressed. In particu¬ 
lar we consider a one-dimensional lattice with a self-affine 
disordered potential, where both localized and extended 
single-particle orbitals were suggested to be present^. 
The second purpose of this paper is to investigate the 
effect of interactions on the Anderson transition. In 
particular, we consider spin-I/2 fermions in the Aubry- 
Andre model with on-site repulsive interactions. This 
model describes the experimental setup recently imple¬ 
mented with ultracold atomic gases in dichromatic opti¬ 
cal lattices by the group of I. Bloch^^. We employ the 
Hartree approximation with temperature-annealed self- 
consistent iterations to determine the phase boundary 
separating the metallic and the insulating ground-states 
in the regime of weak interactions. 

The paper is organized as follows: in Sec. |TT] we intro¬ 
duce the formalism of the MTIS, provide the definition 
of the many-body localisation tensor and describe our 
numerical procedure to compute it. The analysis of the 
one-dimensional Anderson model is reported in Sec. IIIII 
The Aubry-Andre model is analysed in Sec. IIVI and 
the generalised Aubry-Andre model in Sec. |Vl The one¬ 
dimensional Anderson model with long-range correlated 
disorder is studied in Sec. I VII In Sec. IVIIl we address in¬ 
teraction effects in the Aubry-Andre model. In Sec. IVIIII 
we draw the conclusions, focussing on the utility of 
the many-body localisation tensor to identify conductor- 
insulator transitions in disordered systems with and with¬ 
out interactions. 


II. LOCALISATION TENSOR 


where I'k) is the many-body ground state, Rq is the 
a-component of the many-body position operator R = 


Sfci where A is the position operator for particle i, 
with the index i = I,...,V; L is the linear system size. 
We consider ground states of spin-I/2 fermions, with 
N/2 up and N/2 down spins. For noninteracting par¬ 
ticles (or mean-field schemes such as restricted Hartree- 
Focb^S), may be further simplified, givin g^^'^^ = 


det 




, where the matrix 




(with the indices 


j, = 1,2,... N/2) is the overlap matrix whose elements 
are given by 



dr^*(r)e*T’--<^^.,(r), 


( 2 ) 


where (/>j(r) are the single-particle eigenstate spatial 
wave-functions ordered for increasing energies, and r is 
the spatial coordinate. Using this auxiliary quantity, the 
localisation tensor is now defined as^i^ 


\2 ^ ^ 
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where z^^'^ is defined as in Eq. ([T]) with Ra replaced by 
Ra — Rp- In one dimensional systems N = L for half 
filling, and the only component of the localisation tensor 
is the one corresponding to a = /? = x, which is given by 
= -^log|zAr|/27r2. 

In the case of open boundaries, the position operator 
is well defined^i and the localisation tensor may be eval¬ 
uated according to the formula^>2& 

= ((M'|R„R; 3 |vk) - (vI/|R„|vI;)(Mr|R^|vI,))/iV. (4) 


In tune with Kohn’s viewpoint on the origin of the 
insulating behaviour, the authors of Ref. provided a 
quantitative definition of the many-body localisation ten¬ 
sor Xap (the indices a and /3 indicate spatial directions) 
rooted in the modern theory of polarisation. This quan¬ 
tity measures the degree of localisation of the particles 
in the many-body ground state and permits to discrim¬ 
inate metallic and insulating phases. In metals is 
expected to be divergent in the thermodynamic limit, 
whereas it is finite in insulators; thus it defines a many- 
body localisation criterion for the ground state, referred 
to as Kohn’s localisation^. The formula for Xap was orig¬ 
inally obtained from the estimator of the polarization^, 
and can also be derived using a general geometric quan¬ 
tum theory2&. The cases of periodic and open boundary 
conditions need to be treated separately because the po¬ 
sition operator is ill-defined in the former case^i. 

In the case of periodic boundary conditions the locali¬ 
sation tensor is obtained through the auxiliary quantity 
z ^^'^, which for a system of N particles is defined as^i^^ 

( 1 ) 


For a system of independent electrons with N/2 spin-up 
and N/2 spin-down particles, this may be further simpli¬ 
fied to give the squared localisation length 

-^a /3 = J drdr'(r-r')a(r-r')/3|P(r,r')|^, (5) 

where p(r,r') = 2P(r,r') is the one-particle density ma¬ 
trix for a Slater determinant, which is given by2£ 

N/2 

= 2^(/j(r)(/*(r'). (6) 

i=i 

We stress that the length-scale X^x is a many-body locali¬ 
sation length. In particular, it is not simply related to the 
spatial extent of the single-particle eigenstates. For ex¬ 
ample, in the case of noninteracting band insulators, Xxx 
is related to the spread of the maximally localised Wan- 
nier functions^, rather than to the Bloch wave functions. 
Notice that the latter (which are the single-particle eigen¬ 
states) are always extended. There is no simple analogy 
with the Wannier functions for the case of disordered sys¬ 
tems. 

Further insight into the nature of X^^ can be obtained 


= (4'|e*TR=|4-), 
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via substitution of Eq. db]) into Eq. ((S]). In the one¬ 
dimensional case, one obtains the expressioni^: 


2 

^ i=l 

- ( 7 ) 

i=f3 


where x is the single-particle position operator. The first 
sum in this equation is proportional to the second mo¬ 
ment of the single-particle orbitals. The second sum in 
Eq. (O (which in the literature has been refereed to as 
the covariance termi^) originates from the antisymmetry 
of the many-particle wave-function, and would be absent 
in a single-particle analysis. It reflects the correlations 
between different orbitals. These two sums are of the 
same order of magnitude, as confirmed by inspection of 
numerical results. In particular, in the localized phase 
they both contribute to the value of clearly indicat¬ 
ing that this length scale reflects the properties of the 
many-particle wave-function, even in noninteracting sys¬ 
tems. 

In generic insulators, including those with correlations, 
is related to measurable quantities such as the mean- 
square fluctuations of the polarization and the inverse of 
the optical gap via a conductivity sum-rule^. Further¬ 
more, it is related to the spread of the generalized many- 
body Wannier functions, as defined in Ref. 0, which are 
localized in disconnected regions of the high-dimensional 
configuration space, establishing a direct connection with 
Kohn’s theory of the insulating state. 

The formulation to compute the localisation length pro¬ 
posed in Refs. and [^ . and briefly summarized in this 
section, provides a computational procedure to verify 
Kohn’s contention that the many-body ground state con¬ 
tains sufficient information to ascertain whether the sys¬ 
tem is an insulator or a conductor, without recourse 
to the analysis of low-lying excitations. In this pa¬ 
per we provide evidence for a variety of disordered one¬ 
dimensional systems that this is indeed the case. The 
saturation of in the thermodynamic limit is taken to 
signal Kohn’s localisatiort^, whereas its divergence indi¬ 
cates a conductor. 

In our computations we consider both periodic and open 
boundary conditions and employ, respectively, equations 
© and © to compute the localisation length. The 
single-particle spatial wave-functions needed to 

form the many-particle ground-state, are determined 
from full diagonalisation of the Hamiltonian matrix for a 
single spinless fermion using the Armadillo library2£. 



FIG. 1: (Colour online) Half-filled ID Anderson model with 
periodic boundary conditions: Squared localisation length 
X^x (log-scale) as a function of the disorder strength W/t. 
Data for different systems sizes L are shown, and their mu¬ 
tual agreement indicates Kohn’s localisation at all disorder 
strengths. The lines are guides to the eye. Here and in all 
figures the unit of length is the lattice spacing. 

III. ID ANDERSON MODEL 

We consider disordered tight-binding models of nonin¬ 
teracting spin-1/2 fermions defined by the Hamiltonian 

H + ^-C) + VE ^ (8) 

r,cr r,(7 

where r = is the discrete index which labels 

the lattice sites, L is the linear system size, 6 ^, 0 - (^1^) 
is the fermionic annihilation (creation) operator for a 
spin a =t, i particle at site r, and 

corresponding particle number operator. Here and in 
the rest of the article the lattice spacing is used as the 
unit of length, and the (even) total number of fermions 
N is fixed. The hopping amplitudes to the nearest 
neighbours are set by t, Cr is the (random) value of the 
energy at lattice site r, while the parameter W sets the 
strength of the disorder. 

In this section we address, from a many-particle per¬ 
spective, the Anderson model of localisation where the 
on-site energies {cr} are sampled from a uniform proba¬ 
bility distribution in the interval [— 1 ,-|- 1 ]. 

For noninteracting many-particle systems in the ground 
state, the wave function is the Slater determinant formed 
with the lowest-energy occupied single-particle spin or¬ 
bitals. The number of fermions per spin component de¬ 
termines the Fermi energy. In this article, we consider 
the many-particle ground state of N spin-1/2 fermions, 
with N/2 spin-up and N/2 spin-down particles. 

We recall that in one dimensional (ID) systems with 
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uncorrelated disorder all single-particle orbitals are lo¬ 
calised - meaning that they are characterised by an ex¬ 
ponentially decaying envelope - for any nonzero disor¬ 
der strength W. This follows from the one-parameter 
scalingtheory^, and was also shown more rigorously in 
Refs. EH, m According to Anderson’s criterion of local- 
isationii, which is based on the localized shape of the 
single-particle orbitals at the Fermi energy, the system 
should be an insulator at any filling. 

In Fig. [U we show the results for the squared locali¬ 
sation tensor as a function of the disorder strength 
Wjt. The data corresponding to three large (even) lat¬ 
tice sizes with periodic boundary conditions are shown. 
The lattices are half filled, and ensemble averaging of 
the results is performed considering 5—10 realisations of 
the disorder pattern. The localisation length X^x varies 
by a few orders of magnitude as we tune the disorder 
strength. However, it is always finite and system-size in¬ 
dependent, in the whole range of disorder strengths we 
explore, which extends down to the extremely weak dis¬ 
order Wjt = 0.05. These findings constitute a clear sig¬ 
nature of Kohn’s localisation. Also, the variation of A^^, 
with the disorder strength exhibits no sharp features (as 
opposed to the results of the next sections). We veri¬ 
fied that the data obtained using open boundary condi¬ 
tions (not shown) agree with those obtained using peri¬ 
odic boundary conditions. 

Therefore, we conclude that the formalism of the MTIS 
predicts the many-particle ground-state of the ID An¬ 
derson model to be an insulator, in agreement with the 
theory of Anderson localisation and the one-parameter 
scaling theory^. However, in this latter formalism the 
insulating character is attributed to the localised shape of 
the single-particle orbitals in the vicinity of the Fermi en¬ 
ergy, while the MTIS deals with the many-body ground- 
state wave-function. 


IV. AUBRY-ANDRE MODEL 

In this section we consider the one-dimensional Aubry- 
Andre model^^. This is described by the Hamiltonian 
defined in Eq. but with the on-site energies given 
by the incommensurate potential = cos(27rrg -|- 0), 
where g = (v^ -|- l)/2 is the golden ratio, and 6 is an 
(almost) arbitrary phase. This is an archetypal model 
to study Anderson transitions in lower dimensions; it 
has been experimentally realized with ultracold atomic 
gases trapped in bichromatic optical lattices^, and also 
in quasi-periodic photonic lattices^. 

The sinusoidal potential does not display periodicity on 
a finite lattice, and so the Aubry-Andre model is, in this 
sense, disordered. However, this disorder is determinis¬ 
tic, and so is not truly random. In the presence of such 
“deterministic disorder”, as opposed to true disorder, the 
one-parameter scaling theory of Ref. does not apply. 
In fact, it is known that this model hosts a transition 
from a diffusive phase at weak disorder, to a localised 


phase at strong disorder. In the former phase all single¬ 
particle eigenstates are extended over the whole system 
(possibly with the exception of a zero-measure set of non- 
exponentially localised states}^. In the latter phase they 
are all localised if g is a Diophantine number (which is 
the case considered here) and for almost every value of 
9^. We chose 0 = 0. The transition occurs at the critical 
disorder strength Wc/t = 2 . 

In simulations with periodic boundary conditions we 
need to consider system sizes given by Fibonacci num¬ 
bers, so that the potential fits the periodicity of the lat¬ 
tice. The results for the squared localisation length A^x 
of half-filled lattices are displayed in Fig. [2] (top panel). 



0 0.004 0.008 0.012 

1/L 

FIG. 2: (Colour online) Half-filled ID Aubry-Andre model. 
Top panel: squared localisation length A^^, (log-scale) as 
a function of the quasi-disorder strength W/t for periodic 
(PBC) and open (OBC) boundary conditions, and for dif¬ 
ferent Fibonacci lengths L of the chain; the vertical line in¬ 
dicates the critical disorder strength Wc/t = 2, separating 
extended from localised states. Bottom panel: Scaling of 
Xxx with the inverse system size 1/L in the conducting phase 
Wc/t < 2. The continuous curves represent the power-law 
fits X~x = cL ' (c and 7 are fitting parameters). X~x di¬ 
verges with the exponent 7 « 1.14 for PBC (left axis) and 
with 7 = 1.008(5) for OBC (right axis). 
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both for periodic and open boundary conditions. A sharp 
variation of A^x occurs in the close vicinity oiWjt = 2 . 
For W/t > 2, the localisation length is finite and does 
not depend on the system size. This is a signature of 
Kohn’s localisation. Instead, for W/t < 2, a very rapid 
increase of Axx with the system size is observed, pos¬ 
sibly indicating a metallic phase. In order to confirm 
this supposition we perform a detailed analysis of the 
finite-size scaling of Various datasets obtained in 
the regime IF < 2 are shown in Fig. [5] [bottom panel]. A 
best-fits analysis indicates that these data are accurately 
described by the (empirical) power-law fitting functions: 
X~^ = cL~'^, where c and 7 are the fitting parameters. 
The exponents obtained from the fitting procedure are 
7 = 1.135(2), 1.151(2), 1.14(1) (for W/t = 0.2,0.5,1), 
and 7 = 1.008(5) {W/t = 1) for periodic and open bound¬ 
ary conditions, respectively. This fitting function pre¬ 
dicts a divergence of the many-body localisation length 
in the thermodynamic limit, providing a clear indication 
of metallic behaviour. The divergence occurs both for 
periodic and open boundary conditions, but it is more 
rapid in the former case. 

It is worth noticing that in the insulating phase the val¬ 
ues of Axx obtained using periodic and open boundary 
conditions are indistinguishable within our numerical ac¬ 
curacy. This independence from the type of boundary 
conditions is indeed expected for insulators, since in these 
systems the localisation lengths (and the polarisation) 
are bulk properties, as opposed to metals where they de¬ 
pend on the size of the system. 

The analysis of the Aubry-Andre model in the frame¬ 
work of the MTIS provides a clear signature of the metal- 
insulator transition at W/t = 2, in agreement with the 
Anderson criterion of localisation, which also predicts 
a phase transition at the same disorder strength since 
the single-particle orbitals change from extended to lo- 
calised^^. We point out that we also performed a similar 
analysis of the Aubry-Andre model at different lattice 
fillings in the regime 0.1 < N/{2L) < 0.9, without ob¬ 
serving measurable shifts of the critical point. This is 
also expected following Anderson’s criterion of localisa¬ 
tion, since the single-particle spectrum of this model does 
not host mobility edges^^. 


V. GENERALIZED AUBRY-ANDRE MODEL 



W/t 



FIG. 3: (Colour online) Generalized ID Aubry-Andre model 
with next-nearest hopping t 2 /t = 0.5, at half filling. Top 
panel: squared localisation length (log-scale) as a func¬ 
tion of the quasi-disorder strength W/t. Data for periodic 
(PBC) and open boundary conditions (OBC) are shown, for 
different chain lengths L. The gray vertical stripe indicates 
the approximate location of the critical point between the 
metallic and the insulating phases. Bottom panel: scaling of 
inverse squared localisation length l/A^^, with the inverse sys¬ 
tem size l/L at quasi-disorder strength in the metallic phase 
W/t = \ and in the insulating phase W/t = ?> (blue dashed 
lines are a guide to eye). The continuous red curve represents 
the power-law fitting function X^^ = cL with the best-fit 
parameter 7 = 1.19(2). 


Hopping processes beyond nearest-neighbour sites can 
dramatically alter the localisation properties, even caus¬ 
ing the occurrence of single-particle mobility edges when 
none existed in the absence of such effects^^. In this sec¬ 
tion, we consider the generalised Aubry-Andre model, 
including next-nearest neighbour hopping. With this 
modification, one obtains the Hamiltonian H' = H + 

^2 Yl,r,a{H-,a^r+ 2 ,a + where H is defined in Eq. ([ 8 |), 
<2 is the energy associated to hopping to next-nearest 
neighbours, and the on-site energies are defined by the 
same incommensurate sinusoidal potential of the native 


Aubry-Andre model considered in the previous section. 
This generalised Aubry-Andre model was studied in 
Ref. The analysis of the single-particle spectrum 
based on calculations of the inverse participation ra¬ 
tio (which measures the spatial extent of the single¬ 
particle orbitals) presented evidence of the presence of 
mobility edges in a certain regime of disorder strength 
W. The location of the mobility edges was found to 
vary with W. As in previous sections, here we analyse 
the many-particle ground-state of the generalised Aubry- 
Andre Hamiltonian in the framework of the MTIS. We 
consider various lattice filhngs, varying from vanishing 
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FIG. 4: (Colour online) Zero-temperature phase diagram of 
the generalised Aubry-Andre model with t 2 /t = 0.5, as a 
function of filling p = N/(2L) and quasi-disorder strength 
W/t. The black points indicate the phase boundary separat¬ 
ing the metallic phase (left) from the insulating phase (right) 
obtained within the MTIS. Our results may be compared with 
the critical point extracted from the inverse participation data 
of Ref. (blue dashed curve). 


density to full filling. We fix the next-nearest neighbour 
hoppin g a t ± 2/1 = 0.5, a value which was also considered 
in Ref. An illustrative example of the dependence of 
the squared localisation length as a function of W is 
shown in Fig. [3] (top panel), where the datasets corre¬ 
spond to half-filled lattices of different sizes L. Here too, 
as in the case of the native Aubry-Andre model, a sharp 
variation of occurs at a finite disorder strength Wc- 
For W > Wc, A^x finite and system-size independent, 
indicating Kohn’s localisation. Instead, for W < Wc, 
rapidly increases as L increases. In order to assert 
whether in this regime the ground-state is metallic, we 
analyse the finite-size scaling of A^x: Fig. [3] [bottom 

panel]). The scaling of the squared localisation length 
with the system size L turns out to be accurately de¬ 
scribed by the empirical fitting function A~x = cL~"A 
where c and 7 are fitting parameters. At W/t = I, the 
best fit is obtained with 7 = 1.19(2). This scaling be¬ 
haviour clearly indicates a divergence of the localisation 
length, which is a signature of metallic behaviour. 

In order to approximately pinpoint the phase boundary 
between the metal and the insulator, we determine the 
largest disorder strength where Axx clearly diverges in the 
thermodynamic limit, and the smallest value of W where 
it is system size independent, within numerical accuracy. 
This allows us to provide a (rather narrow) interval con¬ 
taining the critical disorder strength Wc- For the case 
of half-filling we obtain Wc/t = 2.275 ± 0.125. This is 
displayed in Fig. [3] (top panel) as a gray vertical stripe. 
By performing a similar analysis for different fillings, we 
obtain the zero-temperature phase diagram as a function 
of disorder strength and filling factor p = N/{2L) (see 


Fig. a. The phase boundary separating the metallic 
and the insulating phases varies rapidly with the filling. 
Interestingly, these variations are non-monotonic: start¬ 
ing from the zero-density limit, Wc first decreases as the 
filling increases, then it rapidly increases when the filling 
is p > 0.5. 

These findings obtained within the MTIS can be com¬ 
pared with the prediction based on the Anderson crite¬ 
rion of localisation. We extract the location of the single¬ 
particle mobility edge from the contour plot data of the 
inverse participation ratio provided in Ref. [ 2 ^ This pro¬ 
cedure is based on the digitalisation of the colour-scale 
shown in Ref. [ 2 ^ and so it entails some approximations. 
Vanishing values of the inverse participation ratio indi¬ 
cate extended single-particle orbitals, while finite values 
indicate localised states. The critical filling factor is ob¬ 
tained when the Fermi energy reaches the mobility edge. 
Notice that in Ref. [H only the single lattice size L = 500 
was considered, without analysing the finite-size scaling 
behaviour. From the scattering of their data for L = 500, 
we estimate the indeterminacy on the extracted critical 
filling to be close to 10%. Therefore, performing a precise 
quantitative comparison between our data and those of 
Ref. [ 2 ^ may not be completely justified. However, we 
see from Fig. 2] that the overall agreement is good. In 
particular, certain important feature of the ground-state 
phase diagram are predicted by both theories. First, in 
the low-filling limit both theories predict the critical dis¬ 
order to be Wc/t ~ 3, significantly larger than in the na¬ 
tive Aubry-Andre model (where Wc/t = 2 ). Second, the 
location of the phase boundary has large non-monotonic 
variations as a function of the filling factor. 


VI. ID ANDERSON MODEL WITH 
CORRELATED DISORDER 

The properties of the Anderson model defined in sec¬ 
tion m are affected in a non-trivial manner by the pres¬ 
ence of spatial correlations in the disorder pattern, in 
particular if the correlations have a long-range charac- 
te ri9i24,36,37 ^ The effect of long-range spatial correlations 
has been investigated in several studies, considering in 
particular disorder patterns characterised by a power-law 
spectral density S{k) oc k~°‘, where S{k) is the Fourier 
transform of the two-point correlation function (e^ £,■'), 
and the brackets (■ • •) indicate spatial averaging. The 
value of the exponent a determines the extent of the 
spatial correlations. The case a = 0 corresponds to un¬ 
correlated disorder. For a > 2 one has energy sequences 
with persistent incremental^. A disorder pattern with 
power-law spectral density can be constructed using the 
following equation^!: 

L/2 

er = ^ (A:““[27r/L]^““)^^^ cos {2'Krk/L -\- 4>k), (9) 

fc=i 
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FIG. 5: (Colour online) Half-filled ID Anderson model with 
long-range correlated disorder and PBC: Scaling of ensemble- 
averaged inverse squared localisation length as a func¬ 

tion of the inverse chain length 1/L (log-log scale), for various 
values of the exponent a characterising the disorder spectral 
density. For a = 3.0, 4.0, 5.0 our simulations indicate a diver¬ 
gence of with system size, which is accurately described 
by the power-law fitting function = cL ' (c and 7 are 
fitting parameters), shown as red solid lines. For a = 0.5, 1 
the results suggest a saturation of A^a, in the thermodynamic 
limit; the blue dotted lines indicate linear fits. 


where (j)k (with k = l,...L/ 2 ) are random phases 
sampled from a uniform distribution in the range [0,27r]. 
In our calculations, we shift the on-site energies in 
order to obtain a disorder pattern with zero mean^l: 
Cave = = 0- Also, in Order to curtail the growth of 

the disorder fluctuations as the system size increases, it 
is necessary to fix the magnitude of the variance of the 
on-site energies &t i^r — eave)^ = Ij by appropriately 
rescaling the on-site energy distributional. Notice 
that, in the notation of Eq. ([5]), the disorder-strength 
parameter is fixed at W/t = 1. This value is not equal to 
the maximum amplitude of the on-site random potential. 
This model of correlated disorder is non-deterministic, 
and so it differs in nature from the deterministic disorder 
of the Aubry-Andre models described in Secs. IIVI 
and 13 However, its properties are also qualitatively 
different compared to the uncorrelated Anderson model 
of section IIIII In fact, a renormalisation group studyll 
predicted that for large exponents a > ac = 2 the 
single-particle orbitals become extended in a finite 
portion of the energy spectrum close to the band centre. 
This is in sharp contrast with the case of uncorrelated 
disorder, where all single-particle orbitals are localised in 
one-dimensional systems^Sr— . It is worth emphasizing 
that the rescaling of the on-site energies was found to be 
crucial for the occurrence of the single-particle extended 
states^^— . 

Fig. [^reports the squared localisation length in half- 
filled chains for various values of the exponent a. The 
data points correspond to ensemble averages obtained 


using from 10 to 200 realisations of the disorder pattern. 
The error bars represent the standard deviation of the 
population (instead of the estimated standard deviation 
of the average) since, as pointed out in Ref. |4l|, in the 
presence of long-range correlations sample-to-sample 
fluctuations survive in the thermodynamic limit. 

As seen from the log-log scale of Fig. [51 the datasets 
corresponding to large exponents a > ac display a clear 
power-law divergence of the localisation length A^x with 
the system size, indicating metallic behaviour in the 
many-particle system. Similarly to previous sections, we 
fit the disorder-averaged data for a > 2 with the fit¬ 
ting function = cL~'^. The best fits are obtained 
with the exponents 7 = 1.26(3), 1.130(4), 1.127(4), for 
a = 3.0,4.0,5.0, respectively. Notice that the latter two 
data are close to the values found for the Aubry-Andre 
model in the metallic phase. In contrast, for lower values 
of a < 2 A saturates with the increase of system size. 
This is confirmed by performing linear fits which extrap¬ 
olate to finite values in the thermodynamic limit. These 
findings are consistent with the single-particle analysis of 
Ref. Indeed, at half filling the Fermi energy is close 
to the band-centre, where the extended single-particle or¬ 
bitals have been predicted to occur. Notice that, accord¬ 
ing to Ref. nil, ensemble-averaging causes the transition 
between single-particle localised and delocalised states to 
morph into a cross-over. However, our aim here is only to 
identify the two many-particle regimes with metallic and 
insulating phases, without focussing on the precise loca¬ 
tion of the transition point. In fact, we have checked that 
the two phases (with, respectively, diverging and saturat¬ 
ing localisation lengths) can be unambiguously identified 
also by analysing the data corresponding to single reali¬ 
sations (not shown). 

Our study, from the many-particle perspective of the 
MTIS, substantiates the claim made in Ref. vis-a-vis 
the occurrence of metallic states in ID. 

VII. INTERACTING AUBRY-ANDRE MODEL 

Understanding the intricate interplay between disorder 
and interactions in many-fermion systems is an outstand¬ 
ing open problem^. In particular, it is still unclear how 
an Anderson localized system is affected when electron- 
electron interactions are include d^^’^^ . The MTIS is 
clearly a promising approach to address this problem, 
given that it allows us to describe both noninteracting 
and interacting insulators within the same formalism. 
Here we investigate the effect of weak repulsive interac¬ 
tions in the Aubry-Andre model. The Hamiltonian we 
consider is defined as 

Hint. = iL + U ^ (10) 

r 

where H is the noninteracting Hamiltonian defined in 
Eq. with the incommensurate disorder pattern 








u/t 


FIG. 6: (Colour online) Interacting Aubry-Andre model at half filling. Left panel: Squared localisation length for a periodic 
chain of length L = 610 for various interaction strengths U ; the jump in signals the conductor-insulator transition. Right 
panel: Paramagnetic ground state phase diagram showing the dependence of critical disorder Wc separating insulating and 
conducting phases as a function of weak interactions U. Inset shows scaling of AjT^ with inverse system length L for U/t = 0.3 
at two fixed disorder strengths W/t = 2.3, 2.7; the scaling behaviour clearly signals an insulator for the latter and a conductor 
for the former. 


employed in Secs. IIVI The parameter 17 > 0 character¬ 
izes the interaction strength, while is the spin- 

up (spin-down) density operator at site r. This Hamil¬ 
tonian describes the experimental setup recently imple¬ 
mented by the group of Bloch—. The experimenters cre¬ 
ated quasi one-dimensional tubes with two standing laser 
waves along the axial direction. One of the two lasers has 
a period which is incommensurate with the other. This 
creates the deterministic disorder pattern er characteriz¬ 
ing the Aubry-Andre model. The interaction strength U 
can be tuned employing a Feshbach resonance. 

Our goal is to determine the zero-temperature phase di¬ 
agram of the Hamiltonian Eq. (fTUl) at half-filling N = L, 
in the regime of weak interactions. We restrict our 
analysis to paramagnetic phases (allowing charge-density 
waves), and we determine the phase boundary separat¬ 
ing the metallic and the insulating ground-states. In the 
regime of relatively strong disorder and weak interactions 
the Hartree approximation is expected to provide reli¬ 
able result a^^d^ . Within this approximation, the Hamil¬ 
tonian is simplified using a mean-held decoupling of 
the interaction term, obtaining: 


Hint. 

I 


i 

+ U 

i 

-U 


( 11 ) 


where H'^ is that part oi H = -|- correspond¬ 

ing to spin (j =t, The densities are obtained via 
a self-consistent iterative procedure based on the equa¬ 


tion (nr,cr) = Qr^a*Qr^aj where is the ma¬ 

trix of eigenvectors of the mean-held Hartree Hamilto¬ 
nians Paramagnetism is enforced by setting 

Special care has to be taken in or¬ 
der to ensure that the iterative procedure converges to 
the true ground-state; following Ref. we implemented 
an annealing scheme where a hctitious temperature pa¬ 
rameter is gradually reduced down to zero temperature. 
This temperature-annealing scheme is combined with the 
standard damping of the density prohle provided by each 
iteration. The Hartree formalism is based on an ansatz 
that the ground state is a Slater determinant, which al¬ 
lows us to compute as described in Section IP. 

Our main results are presented in Fig. |6l The left panel 
shows the squared localisation length as a function of 
disorder strength W for various interaction strengths; the 
sharp jump in A^^, signals the conductor-insulator tran¬ 
sition. In the conducting phase the localisation length is 
cut off by the system size and is independent of the in¬ 
teraction strength. The inset of the right panel of Fig. [5] 
shows the scaling of X~// with inverse system size for 
U/t = 0.3, at two disorder strengths W/t = 2.3,2.7; 
for the latter value we see that, in the thermodynamic 
limit, X'^^ saturates (signalling an insulator) whereas it 
diverges for the former (signalling a conductor). 

An accurate estimate of the transition point is obtained 
by locating the maximum differential of a polynomial 
function which fits the data of A^^, as a function of W. 
With this procedure we determine the zero-temperature 
paramagnetic phase diagram of the weakly interacting 
Aubry-Andre model, shown in the right panel of Fig. [51 
There is a linear increase of the critical disorder sep¬ 
arating the metallic and the insulating (paramagnetic) 
phases. This indicates that repulsive interactions induce 
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delocalisation. It is worth mentioning that a positive drift 
of the critical disorder strength has been obtained also 
in earlier theoretical studies of the Mott-Anderson tran¬ 
sition in higher dimensions based on dynamical mean- 
field theory^. An interaction-induced increase of the 
localisation length was also previously seen in a one¬ 
dimensional Anderson-Hubbard model within a simple 
mean-field treatment valid in the atomic limit — ; how¬ 
ever, no metallic transition was observed in that study. 
In recent experiments a small linear increase of the criti¬ 
cal disorder strength was observed in a three-dimensional 
disordered optical lattice^ . Moreover, the bosonic inter¬ 
acting Aubry-Andre model was implemented in Ref. 1^ . 
where interaction-induced delocalisation was again ob¬ 
served. Bloch’s group implemented the fermionic inter¬ 
acting Aubry-Andre model^^, and determined the critical 
disorder strength where many-body localisation, which is 
a dynamical phase transition taking place at finite energy 
density, occurs. The critical disorder strength was found 
to increase as a function of the interaction strength for 
weak interactions, echoing our findings for the ground 
state. We propose that the ground state metal-insulator 
transition could be observed in their setup by employ¬ 
ing the technique used in Refs, and where an 
effective force is imposed on the atoms either by shift¬ 
ing the harmonic confinement or by applying a magnetic 
field gradient; such experimental results could directly be 
compared with our phase diagram. 

VIII. CONCLUSIONS 

Developing approaches to locate insulating transitions 
in disordered systems, for noninteracting or interacting 
systems, at zero or finite temperatures, is a central 
problem in condensed matter physics and a subject of 
intense researcbi^r— . In this paper we addressed 
the zero temperature aspects of this problem from the 
perspective of the MTIS. Our findings provide evidence 
that the many-body localisation tensor, a bulk quantity 
measuring Kohn’s localisation, provides a clear signature 


of the insulating transition induced by disorder at zero 
temperature. This was first verified in noninteracting 
one-dimensional models with uncorrelated disorder, de¬ 
terministic disorder due to incommensurate potentials, 
and disorder with long-range correlations described by a 
power-law spectral function. 

In particular, it was verified that the ground state of the 
one-dimensional Anderson model is insulating at any 
disorder strength, in agreement with the one-parameter 
scaling theory^ of Anderson localisation. In the cases 
of deterministic and correlated disorder, we found evi¬ 
dence of metal-insulator transitions, in agreement with 
previous studies on the critical disorder strength and on 
the position of the mobility edge based on single-particle 
theories. 

Finally we investigated the conductor-insulator tran¬ 
sition in a disordered interacting system: Using the 
Hartree mean-field analysis within the MTIS, we found 
that weak repulsive interactions induce delocalisation 
in the paramagnetic ground state of the Aubry-Andre 
model at half hlling, leading to an increase of the critical 
disorder strength required for the onset of insulating be¬ 
haviour. These findings could be observed in cold-atoms 
using available experimental techniques ^5,47,49 ^ 

One very appealing feature of the present approach 
is that it permits to identify the insulating phase 
using only ground-state properties. On the other hand 
alternative approaches to identifying the insulating 
state, such as the Kubo formula for dc conductivity, 
require the knowledge of, at the very least, low-lying 
excited states. Within the MTIS, the knowledge of 
the ground-state many-particle wave-function alone 
suffices, a feature which makes it suitable for large scale 
computational approaches such as, say, quantum Monte 
Carlo simulationsiiii^. 


We acknowledge insightful comments from R. Resta on 
the manuscript, help from M. Atambo with computing 
facilities, related collaboration with G. Gebreyesus, and 
interesting discussions with T. Nguyen and J. Goold. 
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